Nonequilibrium dynamics of polymer translocation: a mean-field model 
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We analyse the dynamics of polymer translocation in the strong force regime by recasting the 
problem into solving a differential equation with a moving absorbing boundary. For the total 
translocation time, rtr, our simple mean-field model predicts that rtr ~ (number of monomers)^'^, 
which is in agreement with the exponent found in previous simulation results. Our model also 
predicts intricate dependencies of rtr on the variations of the pulling force and of the temperature. 
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I. INTRODUCTION 

Understanding of polymer adsorption and transloca- 
tion has important technological and biological signifi- 
cance. Besides well-known applications such as adhe- 
sion and coating, adsorption is responsible for facilitat- 
ing breathing in the lungs [Ij] and translocation for the 
mechanism of transporting DNA and RNA across nu- 
clear pores (e.g., see [3]). Although usually viewed as 
two different phenomena, adsorption and translocation 
are in fact very similar as both processes may be seen as 
having one end or both ends of a monomer chain pulled 
to an adsorbing surface. Indeed, simulation results have 
indicated that for a Rouse chain, both the total translo- 
cation time, Ttr, and the total adsorption time Tad, scale 
like M^-^ where M is the initial number of monomers 
in the polymer globule [1, U] . Recently, it has also been 
recognised that many adsorption and translocation phe- 
nomena are in the strong force regime 0, namely, 
bf is at least a few times of ksT where b is the bond 
length between connected monomers and / is the effec- 
tive force exerted on the monomers in adsorption or in 
translocation. This suggests that polymer adsorption or 
translocation is likely to be a nonequilibrium process. 
Coupled with the well-known observations of ageing Q 
and glassy behaviour 0, H, [§] in adsorbed polymer layer, 
the need for a better understanding of the dynamics of 
polymer adsorption and translocation is in order. Here, 
we present a mean-field model that describes the dynam- 
ics of a polymer globule under translocation in the strong 
force regime. Letting the direction of the translocation 
process be in the negative z direction, we denote the num- 
ber density of monomers in each xy-plane along the z 
axis by (p{t, z). The model is mean- field in the sense that 
the fluctuations in the x and y directions are averaged 
over. In the strong force regime, it is expected that the 
polymer will quickly adopt a "stem-flower" type config- 
uration (c.f. Fig. [1]) We thus set up a differential 
equation with the bottom of the "flower" as a moving ab- 
sorbing boundary. The differential equation based model 
can be numerical solved very efficiently. As our model re- 
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FIG. 1: A schematic of the translocation process. The 
bottom-most monomers are cleared first and the resulting de- 
viation in concentration from the equilibrium condition in- 
duces a pressure gradient that drives the monomers at the 
interface, a distance R{t) away from the pore, to the bottom. 




tains all spatial information along the z-axis, it provides 
better characterization of the translocation process. For 
instance, it allows for the determination of the temporal 
evolution of the centre of mass for the portion of the poly- 
mer to be translocated. Our model confirms the scaling 
law: Ttr ~ M^-^, found in simulation results Q. We also 
obtain novel quantitative predictions concerning how Ttr 
would depend on the pulling force and thermal energy. 



II. BACKGROUND 

We are interested in the strong force regime, i.e., the 
effective force involved is more than a few times of ksT /h 
where b is the monomer-monomer bond length. In the 
case of adsorption, this translates to having the adsorp- 
tion energy, e, to be more than a few kgT . In this regime, 
the polymer will be pulled strongly onto the surface and 
the time scale will generally be much quicker than the 
whole polymer relaxation time, Tr — UJ' '^^^'"^ 

M is the number of monomers in the polymer, m is the 
mass of each individual monomer, 7 is the monomer- 
solution collision frequency. The process is thus far from 
thermal equilibrium and local relaxation dominates. Our 
investigation is therefore fundamentally different from 
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much of the earher works on polymer adsorption at equi- 
hbrium pU . [T^ [T3| . NonequiUbrium polymer adsorption 
dynamics has also recently received much attention and 
most studies focused on the scaling for the adsorption 
time, Ta H, w hich for a Rouse chain is found to scale like 
M^-^ [S, llJ]- For further information on adsorption, we 
refer the readers to a recent review by O'Shaughnessy 
and Vavylonis Q . In the case of polymer translocation, 
interest in the physics community is comparatively more 
recent and most early studies have focused on the low 
force regime, in which the relaxation time is shorter than 
the process of translocation. This allows for the use of the 
Fokker-Planck equation description [Tsl . ITg} or the nucle- 
ation theory [13, • The validity of the above formalism 
has been questioned in [l^ as it is argued that the relax- 
ation time and translocation time are of the same order 
of magnitude. The authors further demonstrate the ex- 
istence of anomalous dynamics in translocation through 
simulations and scaling argument. Anomalous dynam- 
ics in forced translocation was also studied in [20| and 
further explored in [2l| with the use of fractional dif- 
fusion equation. More recently, the total translocation 
time for a Rouse chain, rtr, is investigated in Q where 
the authors argue that the pulling force would only af- 
fect a small portion (a "fold" ) of the polymer at a time 
and starting from this assumption, rtr is found by scaling 
argument to scale like M^-^. However, Sakaue argued in 
[22| that the "folds" picture may only be correct when 
bf/ksT > M^/^. In the paper, the author treats the dy- 
namics of translocation as a tension propagation problem 
and by assuming that each blob is at equilibrium locally, 
a differential equation governing the temporal evolution 
of M under translocation is formulated and then solved 
numerically. The approach is very similar in spirit to ours 
although there is one key difference: we treat the ther- 
mal diffusion and applied force separately while Sakaue 
group them together in the form of an effective force: 
/ = fRo/ksT, where Rq is the initial radius of the poly- 
mer globule. In terms of predictions, for a Rouse chain in 
the strong pulling regime, our model and Sakaue's model 
both give rtr ~ M^'^, in agreement with simulation re- 
sults [J]- On the other hand, our model indicates a much 
more complex relationship for rtr's dependencies on / 
and /cbT. In particular, we find that the scaling law: 
Ttr ~ is only true when f /ksT —> oo, and that rtr is 
found to depend non-monotonically on the thermal en- 
ergy. 



III. A DISCUSSION ON SCALING 



For the problem at hand, the dimensionful parameters 
are: /, /cbT, &, 7 and m, with M as the only dimension- 
less parameter. If we let the total translocation time, 
Ttr, be given by the function: 4>{fTkBT^b,^, 111,^1). We 
can invoke the intuitive Il-theorem '23l| to transform the 



FIG. 2: Temporal evolution of <j}(t,z) with my = 1, 6 = 1, 
/ — 10, D = 1 and AI = 500. The times are in units of tq. 
Notice the gradual lengthening of the stem attached to the 
pore at jz = 0. 




z {in units of b) 



functional dependency into the following form 
bf keT 



Ttr = To$ 



where Tq = ^^^^rp is the single monomer diffusion time 
[ll[ , and $ is now dimensionless and as such, it can only 
depend on the first two and the last parameters, i.e., 

/_&/ ksT 



ksT' mb'^y^ 



,b,y,m,M 



(1) 



Tad = To^> 



.M 



(2) 



The above equation is exact except for, of course, the 
fact that we do not know what $ is. 

In Hi, the author's ansatz for the form of $ is: 



knT 



yksT mb'^y'^ 



,n] ^ const. x'^N^/^ . (3) 



Namely, it is assumed that the second argument in <I> 
is redundant. There is no physical reasoning for this 
particular ansatz. Indeed, we find in this paper that all 
three arguments affect $ independently even in the range 
10 < bf/ksT < 500, which is not described in 



IV. A MEAN-FIELD MODEL 

Before we move on to describing our model, we re- 
mark that in this paper, we ignore all interactions from 
monomers that have passed through the pore. Namely, 
we assume that the pulling force is strong enough that 
the chain on the other side of the pore does not have 
any effect on the chain b eing translocated. The same as- 
sumption is made in 0, [23 | and it renders our analysis 
more similar to the phenomenon of adsorption and so the 
findings here may be relevant to both phenomena. 

The basic observation concerning translocation is that 
as the first monomer of a polymer globule is being 
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FIG. 3: Some properties of (j>{t, z) with parameters as defined FIG. 4: Adorption times with respect M. It is found that 

in the caption of Fig. [J] (a) Center of mass for the portion of Tad ~ M^'^. 
polymer awaiting translocation (in units of h). (b) Number of 

monomers, M{t). (c) The length of stem (in unit of b). , kj.i,m = i,Y.i 
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dragged through the pore (c.f. Fig. [T]), the monomer 
number density at the bottom of the coil gradually de- 
creases. This imbalance in monomer concentration from 
the equilibrium condition will induce an osmotic pres- 
sure that pushes the globule towards the surface. As a 
result, it accelerates the adorption process. Starting with 
this observation, we let ^(t, z) be the expected number 
of monomers at time t and position z. In other words, 
(f)(t, z) has dimension length"-'^. Denoting the total num- 
ber of monomers by M and the center of mass by z, 
we write 0m, j as the monomers distribution at equilib- 
rium. Given any other distribution, 0, we assume that 
the osmotic pressure to be proportional to the difference 
between the current distribution and the distribution at 
equilibrium: (j) ~ 4'M,z (c.f. Ch. VII in [HI). In other 
words, if we ignore adsorption for the time being, the 
temporal equation on the distribution is: 



'dt 



(4) 



where D is the diffusion constant and is assumed to be 
kBT/rw-f. We note that as a deterministic model, the 
above equation does not model diffusion of the whole 
molecules, i.e., z does not vary and as such, the model is 
meant to present the dynamical behaviour at short time 
in comparison to the whole globule relaxation time, = 
M^tq. This assumption is consistent with the parameter 
set we study here as Ttr is always less than 10 percents 

of Ti-. 

We now incorporate translocation into the model. If 
the monomers are not connected, the adsorption process 
may be modelled as a fixed absorbing boundary in the dif- 
fusion equation. But since the monomers are connected 
and as such the pulling force can propagate through the 
chain, the differential equation above becomes a mov- 
ing boundary problem and we have the following model 
equation: 



/ 



dt 



jnjbR 



for z ^ R{t) 
for z > R{t). 



(5) 



where R{t) = max[z : (j){t, z) < b and < z < z{t)], and 
/ is the pulling force. We also maintain that (j){z,t) = 1, 
for < z < R{t), which represents the stem connect- 
ing the pore and the flower (c.f. Fig. [T]). In the above 
equation, R{t) is the moving absorbing boundary with a 
constant absorbing rate — ^^f^j^ ■ The rate equation can 
be obtained from the force-velocity equation: 



bf 

rwyR 



dt 



(6) 



By dimensional analysis, we know that 4>m.z{z) = 
(f>M',z{{M' /MYz) where ly is 3/5 for a chain in good sol- 
vent and it is 1/2 for a chain in 9 solvent [ll|. In other 



words, if we know 



'M(0),0: 



we can obtain all the other 



distributions 4)m,z by simple translation and dilation. 

In summary, we have constructed a differential equa- 
tion model that depends purely on a set of constant 
parameters: /, m, 6, 7, ksT, and a static distribution 
(t>M{o),Q that can be determined once and for all. Eq. [5]is 
the main result of this paper and it can be numerically 
solved efficiently (c.f. Appendices A and B for simulation 
details). We will now focus on the various predictions 
made by our model on the Rouse chain. 



V. RESULTS AND DISCUSSION 

Since the pulling force is strong, the diffusion process 
is not rapid enough to replenish the supply of monomer 
near the pore and so a stem forms naturally as shown 
in Fig. [21 This renders an initial increase in separation 
between the center of mass and the pore as the bottom 
monomers are quickly translocated (c.f. Fig. [3] (a)). To 
understand the scaling law: Ttr ~ AI^-^ (c.f. Fig.[3]), one 
can imagine the scenario where ksT — > 0. In this sit- 
uation, the globule is completely frozen throughout the 
translocation process. Now, since the initial size of the 
polymer globule is ~ by/Mo where AIq denotes the initial 
number of monomers, the stem will be of the same order 
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of magnitude in length for most of the translocation pro- 
cess. Therefore, as a first approximation, one can set up 
the following differential equation as in Eq. [6l 



/ 



dM{t) 



(7) 



and the scaling law will then follow immediately. In 
other words, the scaling is purely due to the fact that 
the monomers being dragged to the pore are on average 
a distance of ^ b^/Mo away. This is a much simpler ex- 
planation of the scaling law than that presented in Q 
and it highlights that the "fold" picture may not be nec- 
essary in explaining the scaling behaviour seen in single 
chain translocation f24j. 

If the pulling force is increased, it is natural to expect 
that Ttr will decrease. This is indeed the case, but devi- 
ation from the expected scaling law: Ttr ^ can be 
seen even for the range 10 < bf /ksT < 500 (c.f. Fig. [3 . 
This is different from the expectation described in 0, 124I . 
In fact, our results suggests that the above scaling only 
holds at the limit //fcsT — > 00 and as such, highlight the 
important role of the thermal energy. 

If the temperature is increased, the diffusion process 
(indicated by the red arrows in Fig.[T]) induced by the os- 
motic pressure is facilitated and one would expect an de- 
crease in Tti-. Although this is generally the case, it is sur- 
prising to see the opposite trend at the low-temperature- 
high- force regime (c.f. Fig (5]). This counter intuitive 
feature may be understood by the fact that at low tem- 
perature, as the force is becomes large, the center of mass 
of the remaining polymer is driven away from the pore 
quickly (as shown in Fig. [3]) and this escape process is 
aided by a slight increase in diffusion as the temperature 
is increased. 

In conclusion, we have formulated a simple mean-field 
model for polymer translocation that captures the effect 
of local diffusion. Our model is capable of confirming the 
scaling law: th ^ M^-^, as observed in simulations [1], 
and predicts an intricate Tti 's dependencies on the pulling 
force and the thermal energy. Our work thus signals an 
interesting new territory that awaits exploration. 
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FIG. 5: Adorption time vs. pulling force. The results indicate 
that the scaling law: rtr ~ is only true asymptotically as 
ksT — > 0, and deviation from it can be observed even for 
ksT = 0.04 and for the range 10 < / < 500. 




FIG. 6: Adorption time vs. thermal energy. Notice that as / 
increases, rtr can become non-monotonic with respect to ksT 
as shown by the peak indicated by the black arrow. 
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Thermal energy, k^T 



where P{x) is 
- 0.6x1° -I- 4a;^- 9.48a;® 



11.21a;^-11.94.x2 + 1.27 (A2) 



for l^l < 1.44 and (t>M,o{x) = otherwise (c.f. Fig. [7]). 
This functional form is used in our numerical integra- 
tion although this approximation step by an closed form 
function is not necessary. Instead, one can formulate a 
lookup table for V^0m,2 from sampling. 



APPENDIX B: 



DETAILS ON SIMULATION 
METHOD 



APPENDIX A: FITTING FOR (f>M,2 

For a Rouse chain in three dimensions with & = 1, wc 
find that the following function is a good approximation 
for (h: 



4>M,ziz) 



M 



■ exp 



P 



(Al) 



In solving the differential equation Eq. [5l we employ 
the finite element method. Namely, we denote (j){tn,Zj) 
by 0" where i„ and Zj are the grid points on time and on 
position. In our simulations. At = 0.0017^^ and Az — 
0.056. Specifically, our algorithm is as follows: 

1. Given M and a time grid and a position grid with 
spacing At and Az, set n = 0, z = 0, and for all j. 
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FIG. 7; Fitting by the expression shown in Eq. Al. 
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set 0" = 4>M,o{zj) where 4>m,o is given in Appendix 
A. Let po = max[j : (jp^ <1 and Zj < z], set s — Zp^, 
R — Az and v = . 



2. For j > p, set f/ij"*"^ as 

3. For po<j< P, set to 1. 

4. Set M as M - vAt and z = J^j ^j</'"^V 
Renormalise '/'j^^ -^^ ^7 re-scahng (t>"^^ ■ 

5. Let p = max[j : 0"^^ < 1 and < z], set i? = 
Zp- s + Az and u = 

6. If M < 1, stop; otherwise, increment n by 1 and go 
back to 2. 
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